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We overview the general status of higher order corrections to Bhabha 
scattering and review recent progress in the determination of the two-loop 
virtual corrections. Quite recently, they were derived from combining a 
massless calculation and contributions with electron sub-loops. For a mas- 
sive calculation, the self-energy and vertex master integrals are known, 
while most of the two-loop boxes are not. We demonstrate with an ex- 
ample that a study of systems of differential equations, combined with 
Mellin-Barnes representations for single masters, might open a way for 
their systematic calculation. 
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1. Introduction 

Since Bhabha's article on the reaction 

e + e" -> e + e" (1.1) 

numerous studies of it were performed, triggered by better and better ex- 
perimentation. Bhabha scattering is of interest for several reasons. It is a 
classical reaction with a clear experimental signature. 

Small angle Bhabha scattering has a specific sensitivity to new physics 
due to its forward peaking structure, and is a prominent luminosity monitor 
at high energies. At LEP, the cross-section predictions were needed with 
an accuracy of the order of about 10~ 3 2 . At the planned linear e + e~ 
collider ILC (see the documents prepared for TESLA (Hj), another order 
of magnitude (10 -4 ) might be needed [1] and realized jSJ. This accuracy 
is in reach of the most advanced Bhabha Monte Carlo programs BHLUMI 
0IZj, NLLBHA [HI El, SAMBHA [TO]. The programs contain higher order 
photonic corrections, both exponentiated soft photon corrections and log- 
arithmically enhanced fixed order corrections (including 2nd order). Con- 
cerning the virtual corrections, the logarithmically non-enhanced two-loop 
contributions (and also radiative one-loop contributions) were not included 
until recently. 

At small energies, large angle Bhabha scattering may be used for lu- 
minosity determination. The KLOE collaboration, e.g., aims at an experi- 
mental accuracy of 0.3 %, compared to a present accuracy estimate of the 
Monte Carlo programs (BABAYAGA QUITO], BHAGENF [TO], BHWIDE 
Q31ITO], MCGPJ [TO]) of about 0.5 % [T7j. Here, the inclusion of two-loop 
boxes is not needed. 

Many calculations of weak corrections have been performed since the 
first complete one- loop calculation |18j . and the technical accuracy is ex- 
tremely high nowadays. A recent one, based on DIANA |19l 1201 121j and 
altalc |22| . demonstrated an agreement of more than 10 digits with another 
calculation [20] . Leading higher order weak corrections due to the top quark 
were included in |23| . where two calculations were shown to agree within 
few per mille. It might also well come out that this accuracy (for the pure 
weak part) is sufficient for applications in the foreseeable future. But, it 
might well be that a two-loop calculation of all the dominant terms of the 
second order weak corrections is of interest. 

The last few years brought considerable progress in the pure photonic, 
virtual two-loop corrections to Bhabha scattering. They were determined 
for the massless case some time ago |24j . and quite recently the under- 
standing was developed on how to transform this into the on mass shell 
renormalization scheme |251 126j . This is an absolutely needed part of the 
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calculation in this approach since the Monte Carlo programs in use rely 
on a finite electron mass, thus regulating the collinear singularities. For 
this, several infra-red properties of the amplitudes had to exploited: ex- 
ponentiation of the IR logarithms; factorization of the collinear logarithms 
into external legs; non-renormalization of the IR exponents. The last of 
the three properties is not fulfilled for diagrams with an internal fermion 
loop. For this reason, the massless calculation has to be combined with an 
explicit calculation of them |27l 1281 129j . In [3U] , the massive calculation of 
the Italian/Freiburg group is summarized. So far, it is gauge invariant but 
not complete since the massive box diagrams are an inevitable part of that 
program. Further, the present combination of the two existing calculations 
(Penin's and Italian/Freiburg approach) seems not yet to match exactly 
what is a physical quantity (Bhabha scattering with several flavors); see 
also |3T] . 

In view of the importance of the problem and the principal interest it 
is quite obvious that a complete massive calculation is to be done yet. A 
substantial last missing part are the massive two-loop boxes. 

2. Master integrals 

The two-loop vertex master integrals have been completely determined 
and published by two groups in 2003 [HH ESI an d 2004 [31] . For details and 
references for the self-energies see [SHI- The masters are available, from both 
groups, from the internet as Form files or MATHEMATICA files [361 I37| . 

Concerning the double boxes, the situation is less developed. A complete 
list of masters is known [331301) but the evaluation is, as was expected, quite 
difficult. Not only that they depend on a more complicated kinematics 
compared to vertices or self-energies, namely on both s and t. Additionally 
one is faced with the fact that up to six masters (for B613m3 1 ) may form 
systems of differential equations, if one wants to go this way. So it may well 
happen that the use of Mellin-Barnes representations for the determination 
of single masters (first successfully used for massless masters [23 EH] ) has 
to be preferred here. In fact, the heroic efforts of V. Smirnov |401 141) give 
strong hints to that. The two-loop Bhabha box diagrams contain Feynman 
integrals with seven propagators (lines) and with numerators, and their 
reduction to a smaller number of master integrals leads to the inclusion of 
masters with a smaller number of lines, and also of masters with dots (with 
lines raised to a power biger than 1). The present knowledge of box masters 
may be summarized shortly: 

1 B613m3 is a master topology with six lines, three of them massive; for details of our 
notations, see 1351 136| . 
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Fig. 1. Two-loop Mis with five internal lines 



• Seven lines: For topology B714ml (first planar) exist two masters, 
and B714ml and B714mlN are known [101111]. For topology B714m2 
(second planar) we need four masters, and the basic master is known 
(with use of a two- fold integral in the constant term) Finally, 
for the three masters of topology B714m3 (non-planar), only the most 
divergent part (~ 1/e 2 ) of the basic master is determined analytically 



• Six lines: The ten masters are not published. 

• Five lines: Analytical results for singularities are known, see next 



The two-loop box masters depend on two variables s and t, or equiva- 
lently on related scaling variables x and y (the mass m = 1). There are five 
master topologies with five internal lines, see Figure ^ 

As experience shows, some of them are already examples of quite compli- 
cated systems of master integrals (Mis), although two cases, namely B512ml 
(one MI [MIES]) and B514m (two Mis [23 El EE] ) are relatively simple. 
B512m2 is a system of three Mis, B512m3 is one of four Mis, while B513m is 
actually a system of even five Mis (see Fig- El) ■ In |35j . analytical results for 
the singularities of topologies B512m2 and B512m3 have been given. Here 
we complete the determination of singularities of the Mis with five lines and 
consider the most complex case B513m. 
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B513md2' B513md3 

Fig. 2. Five Mis for topology B513m. Masters B513md2' and B513md2 are symmet- 
ric, B512m3d2'[s(x),t(y)]=B513md2[t(y),s(x)]. 



3.1. B5l3m: Synergetic use of Mellin- Barnes representations and 

differential equations 

If we try to calculate the masters for topology B513m with differential 
equations (DEs), we meet a system of five coupled differential equations. In 
fact, one of the dotted masters (B512m3d2) appears twice (with interchanged 
arguments). 

Fortunately, we can use another method, based on Mellin-Barnes pre- 
sentations, to solve some of them separately. We found it most efficient to 
solve the case B513md2, where one eliminates two masters at once. Details 
of the calculation will be given elsewhere |42j. The result is: 

B513md2 = X 9 . H[0,x] 

e A 4(1 — x z ) 



+ 



1 x 



2{l + y 2 )H[0,x)H[0,y] 



e 4(1 - x 2 )(l - y 2 ) 
- (-1 + y 2 )(6 + 6H[-1, 0, x] - 4H[0, 0, x] - 2H[1, 0, x]) 
+ 0(1), (3.1) 

Then, up to (and including) order 1/e, we are left with a system of only 
three DEs. The differential equation for B513md3[-1] is: 

dB5l3md3[-l] 1 + x 2 yH[0,y] 

Wx = x0^) mi3md3[ - 1] ~ (I - x 2 )(l - y 2 Y (3 - 2) 

with s = —(1 — x) 2 /x, t = —(1 — y) 2 /y, and can also be solved separately. 
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The solution is: 

B5l3md3[-1] = — * V H[0,x}H[0,y}. (3.3) 

(-1 + ar)(— 1 + y l ) 

For our notations in terms of HPLs and generalized HPLs see |43| and 
references therein. 

Finally we have to determine the singularities of B513mdl. They can be 
found from algebraic manipulations (see [II])- This relies on the fact that 
the basic integral, B513m, and the integrals with numerators are nonsingu- 
lar. In order to use this, we replaced in the defining equations for B513mdl 
the dotted diagrams by those with numerators. The solution is: 



B513mdl 



1 1 



e 2 8x(-l + y)(l + y) 3 

\-l + x) 2 y(-l + y 2 + 2yH%y}) 

1 y 



e24z(l + x)(-l + y)(l+y) 3 ) 
x [6(-l + x - x 2 + x 3 )H[0, x]{-\ + y 2 + 2yH[0, y}) 

- 6(1 + x)(-2 - 2x 2 + 2y 2 + 2x 2 y 2 + y£ 2 - 2xy& 
+ x 2 yt,2 + 2(-2x — y + 2xy - x 2 y - 2xy 2 
+ (-1 + x) 2 yH[-l, -y] + 3(-l + x) 2 yH[-l, y])H[0, y] 

- 6(-l + xfyH [0, -1, y] - 4yH[0, 0, y) + 8xyH[0, 0, y] 

~ 4x 2 yH[0, 0, y] + 2yH[0, l,y] - 4xyH[0, 1, y] + 2x 2 yH [0, 1, y\) 
+ O(l), (3.4) 

We can easily check that the solutions (J3.4D and ()3.1I) fulfill the original dif- 
ferential equations (not reproduced here for space reasons). The knowledge 
of the solution ()3.1|) . obtained by the Mellin-Barnes method, has been used. 
So, our combined approach is also a cross check for the analytical result 
obtained with the Mellin-Barnes method. Finally, we checked that (|3.4|) . 
()3.3|) . and (|3,1|) are also in agreement with a numerical approach based on 
sector decomposition calculations. 

4. Summary 

We gave an overview on status and recent progress for the higher order 
corrections for Bhabha scattering. The small angle Monte Carlo programs 
are approaching the accuracy standard of 10 -4 once this is needed. 

The results of the determination of the pure constant (non- logarithmically 
enhanced) virtual two-loop corrections seem to stay perhaps negligibly small 
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(at this accuracy level and after isolating from them not only the mass-, but 
also the kinematically enhanced terms) in the kinematical region of interest, 
but this is not finally numerically fixed. 

In order to have a complete massive calculation (though in the limit of 
small electron mass), one has to determine yet the bulk of two- loop boxes. 
This is underway. 
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